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Abstract 

2-point correlators of vector and axial currents, obtained 
from r-decay data, are studied within the framework of 
perturbative QCD and Operator Product Expansion. Var- 
ious sum rules, obtained from Borel transformation of the 
correlators in complex plane, are used to separate the con- 
tributions of different operators from each other. The anal- 
ysis confirms the Q 2 -dependence of the correlators in the 
space-like region, predicted by QCD+OPE. However the 
operator values are found to be in certain disagreement with 
the estimations, obtained from other data. 

1 OBJECTIVES 

Precise measurements of vector V and axial A spectral 
functions in hadronic r-decays by ALEPH [^J and OPAL 
|Q| collaborations provide us with possibility to test vari- 
ous QCD aspects. Perturbation Theory (PT) and Opera- 
tor Product Expansion (OPE) are the most well-established 
ones. Here we shall compare theoretical predictions with 
the data within the framework of sum rules. Particular de- 
tails of this analysis can be found in 

The 2-point correlators of charged vector and axial- 
vector currents 



J =V 1 A: V^ = wy^d , A^ = uj^d 
can be parametrized by 2 polarization functions H(q 2 ) 

i J e iqx (j fl {x)J v (0)^)dx = 



(1) 



For q 2 — s > they have imaginary parts, the so-called 
spectral functions 



vi/ffli(s) = 27rImn^L(s + iO) 



(2) 



which have been measured from hadronic r-decays for < 
s < m 2 , the plots can be found in 0]. The scalar axial 
polarization function Iiy is basically saturated by single 
pion decay channel. Its imaginary part o,q is delta-function, 
which can be easily separated from ai jv\ . 

It turns out to be convenient to consider the sum and dif- 
ference v\ ± a\ instead of v\ and a\ separately. Indeed, the 
sum v\ + ai is known with better accuracy, while the dif- 
ference v\ — ai does not contain perturbative terms in the 



massless quark limit. The QCD expressions for appropriate 
polarization functions can be written in the following form: 
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The 2fc-dimensional constants are the vacuum expec- 
tation values of the operators, constructed from the quark 
and gluon fields [^J. They have been computed up to di- 
mension D = 8. The numerical values of cannot be 
determined within the perturbation theory. 

Obviously the expressions are not valid for all val- 
ues of s. Exact polarization operator H(q 2 ) is known to be 
an analytical function of s = q 2 with a cut along positive 
real semiaxes. So it is convenient to study the QCD predic- 
tions in the whole complex s-plane. These series are 
not valid for small |s|, where effective degrees of freedom 
are hadrons rather than quarks. Moreover, the higher loop 
perturbative terms in (Q) have unphysical cut starting from 
some s = — Qq < 0. The OPE series with finite number of 
operators does not have a cut along positive real semiaxis, 
but has very singular behavior at s = 0. Based upon these 
speculations one may draw schematic Figure [j], displaying 
the region of validity of the series (U[|). 

Another drawback of QCD is that these series are likely 
to be asymptotic, i.e. divergent for any fixed s. The way to 
deal with divergent series is to apply Borel transformation 



Bm 2 II = pert, terms + 
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Figure 1 : Region of validity of perturbation theory and op- 
erator product expansion 



which improves the convergence by suppressing the higher 
terms. It is not clear, whether it improves the perturbative 
series, which is an expansion in inverse powers of In (— s), 
rather than s itself. However the expansion in 1/ In M 2 
might be convergent. 

The primary goal of the investigation is to find the nu- 
merical values of input theoretical parameters, such as 
a s (m 2 ) and few operators of lowest dimensions. We shall 
compare the QCD result for the Borel transform (Q) of the 
series with the experimental values, computed by ex- 
ploiting the analytical properties of exact polarization func- 
tions. In order to separate the operators from each other, we 
shall consider the Borel transform ^ at complex values 
of the argument M 2 e 1 ^. This can be alternatively under- 
stood as the Borel transformation applied to the polariza- 
tion function, taken at the angle tj> w.r.t. the real negative 
semiaxes in the s-plane, see Fig. [I] We shall also try to 
find the lowest value of the Borel mass M 2 , at which the 
comparison of QCD to the experiment can be made. 

2 V-A SUM RULES 

We start the analysis from the V — A case which is 
purely nonperturbative. The dispersion relation for the dif- 
ference of polarization functions does not need subtractions 
and is written in the following way: 



n«( s )-n^( s ) 



(!)/ 



(6) 



The last term is the kinematic pole which is specific feature 
of axial currents. Indeed, the r.h.s. has appropriate asymp- 
totics: at s — ► it matches the chiral theory prediction, 
while the expansion at s — > oo starts from the operator of 
dimension D = 4, as it should be. 

Applying the Borel transformation to (^), one gets the 
following sum rule: 
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One may estimate the numerical values of the operators 
O v ~ A up to dimension 8 from other data: 



0\' A = 2(m u + m d ){qq) = -f 2 ml 

2 



o, 



V-A 



negligible at s ~ 1 GeV 

64 
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y 



87ra s mo {qqy « 2 x 10~ 3 GeV°, (8) 



where 



qGq 



0.8±0.2GeV^ 



(9) 



has been found from barionic sum rules In the nu- 
merical estimation we assumed m u + nid = 12 MeV and 
a s = 0.5 at 1 GeV 2 . The factorization hypothesis was used 
in order to bring the operators Oq,s to the form (^). It has 
internal theoretical ambiguity ~ 1/N 2 among the D = 8 
operators, see [S], 

QCD corrections to the operators Oq have been com- 
puted in [Q. They turn out to be large and may increase the 
effective contribution of the D = 6 operator by about 50%: 



O, 



V-A 



— 9~ ' ' 
-3 x 10~ 3 GeV 6 



^f-ln^ 



C6 



(10) 



The coefficient cq is ambiguous: two essentially different 
choices were presented in [Q]. In the numerical estimation 
( |l(i| ) we used more moderate one cq — 89/48. 

One sees, that the r.h.s. of (Q) has leading term f 2 and 
relatively small (but interesting) contributions of 0$.$ at 
M 2 > 0.5 GeV 2 . One way to kill f 2 is to differentiate (@) 
by M 2 . This however inevitably increases the errors of the 
experimental integral. It seems more effective to perform 
another trick: one substitutes complex Borel mass M 2 e 1 ^ 
into (R) and takes imaginary part of it. The result is: 
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Let us consider the angle (f> = tt/3. The operator Og 
disappears from the r.h.s. of (jllj) and only Oq is important 
in this case. The l.h.s. of ( fi"l| ) is shown in Fig. ||a as shaded 
area (the upper integration limit is m 2 ., since there are no 
data beyond this point). The errors have local minimum at 
the point M 2 = 0.8 GeV 2 . It happens because the sin (...) 
in the integral has zero at s = m 2 and thereby suppresses 
large experimental errors. At this point we determine the 
operator A and plot the r.h.s. of ( pT| ) with this value in 
Fig. ^a as dash line. 

Second interesting angle is cf> — it/ 4. Both 6 and O s 
contribute, but the next term with O 10 disappears. This 
means that one may go to lower values of M 2 in order 
to reduce the experimental uncertainty. Indeed, as can be 
seen from Fig. ^Jp, the agreement can be achieved down to 
M 2 = 0.4 GeV" in this case. At this point we obtain the 
most accurate value of the operator 0$- 

The result of the fit: 



V-A 
6 

V-A 



-(6.8 ±2.1) x 10~ 3 GeV 6 
(7 ± 4) x 10~ 3 GeV 8 



(details of the fit and error estimations are discussed in [pQ) . 
The result for Oq is twice larger than our estimation (|l0|). 
It might have different explanations: overestimated rn u + 
md, failure of factorization, large a s corrections are the few 
ones among them. But the mass to 2 ,, obtained from ( |l2| ) is 
in agreement with ([)]). 
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Figure 2: Sum rule (Ml) for <fi — ir/3 (a) and <fi — n/4 (b). Dash lines display OPE prediction with operators (113). 



3 PERTURBATIVE SERIES 

Before analyzing V+A sum rules, we outline the basic fea- 
tures of perturbative series. The QCD coupling a = o. s /tt 
is a function of the scale Q 2 , determined by the renormal- 
ization group equation: 



da 



dhxQ' 



k>0 



(13) 



where the factors (3^ have been computed up to 4 loops in 
MS scheme [g|. In particular /? = 4/9, = 4, fa = 
10.06 and /3 3 = 47.23 for 3 flavors. The solution of RG 
equation is 



In 



a(p 2 ) /^( a ) 



(14) 



Since the integral is convergent at oo for any fixed order (at 
least with positive (3k), the coupling a(s) has unphysical 
singularity at some negative s = —Qq, see Fig. [j]. The 
properties of the solution of RG equation can be understood 
by viewing on Fig. ^|. 

The polarization function is obtained by integrating the 
Adler function, which is finite and has been computed up 
to N 3 LO term in MS [01: 
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where K2 = 1.64 and K3 = 6.37 for 3 flavors. In our 
calculations we shall take the theoretical uncertainty equal 
to the contribution of the last term in (|l5|), ztK^a 3 . Since 
we do not know K4, it is reasonable to use only 3-loop 
approximation also for the /3-function in ( 13 ). 

The polarization function constructed in this way has un- 
physical cut from s = — Q 2 to s = 0. It is an obvious indi- 
cation of QCD inapplicability at low \ s\. However there are 
certain attempts to construct the perturbative functions with 



appropriate analytical properties on the whole s-plane, for 
instance by constructing an analytical QCD coupling with 
help of dispersion relation [|Tot (subtractions assumed): 



a s (s) a 



Ima s (s') , 
■ ds 



A 2 



O VM-s/A 2 ) A 2 + sJ ' (16) 

This way is not unique: one may write down the same 
dispersion relation for the polarization function II(s) an as 
well. At the NLO level the result will be the same as the 
substitution of ( |Toj ) into (|i~5|), but it is not the case for higher 
terms. In general, the purely logarithmic terms in analytic 
approach are the same as in conventional QCD, but the 
power terms are different, and there appears D — 2-like 
term ~ 1/s, absent in canonical OPE. 

QCD must give correct value of the hadronic r-decay 
branching ration R T ~ 3(1 + ), which is measured with 
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Figure 3: Real and imaginary parts of a^|g(s)/7r as exact 
numerical solution of RG equation ( pj| ) on real axes for 
different number of loops. The initial condition is chosen 
a s = 0.355 at s = —m 2 . Vertical dotted lines display the 
position of the unphysical singularity at s = — Qq for each 
approximation (4 — > 1 from left to right). 
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Figure 4: Correction 8^ versus a s (m 2 ) and a s {m 2 z ) in 
conventional and analytic approach in 3-loop approxima- 
tion. 



rather high accuracy. It is also weakly sensitive to the non- 
perturbative power corrections. The perturbative fractional 
correction is given by well-known formula (e.g. KW): 



1 + (5 (0) = 1.206 ±0.010 
= 2ni <f *L ( 1 - -L 

J | s |=m? 
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(17) 



Notice, that the circle integral includes the contribution 
of unphysical cut, while in any analytical approach it is 
thrown away. The numerical results for ( ]17[ ) in both ap- 
proaches are shown in Fig. || It is seen, that analytical 
scheme predicts very large a s (m 2 z ) (at Q 2 = m 2 z the dif- 
ference between both approaches is not important) and, 
therefore, fails to agree with other data. So we shall not 
consider it anymore. 
It follows from Fig. ^] 



a s {m 2 T ) = 0.355 ±0.025 



(18) 



The error includes the theoretical uncertainty iK^a 3 of 
Adler function. This result is 3-loop, the 4-loop result with 
the estimation K4 = 25 — 50 would give us slightly less 
value within the error range. We notice also, that there 
is second point on Fig. |] where the conventional curve 
crosses the experimental band. However it is unstable un- 
der changes of various perturbative input parameters and 
prescriptions and cannot be considered as reliable one. 



The operators O v ^ 
sates. In particular. 
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include purely gluonic conden- 
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The D = 4 gluonic condensate has been found from char- 
monium sum rules 
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The D = 6 operator contains gluonic condensate ~ (G 3 ), 
which is not known. The quark contribution after factor- 
ization get the form: 



O, 
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(1.3±0.5)xlO" 3 GeV 6 (20) 



For numerical estimation we used our V - A fit @ and 
added additional error which might occur due to incom- 
plete cancellation of two relatively large term in the sum 
V + A after factorization. The D = 8 operator cannot be 
obtained from other data, but we estimate its upper limit as 
Details given in [ 



\Ol +A \ < IO- 3 GeV 8 



are essentially smaller than 0^$ A and perturbative terms 
dominate here. 

Now let us define the Borel transform of the polarization 
function Wy+A'- 
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The perturbative part B pt is computed numerically: 

ds 



B p t(M ) = in 
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Hpt(s) 
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The integration contour goes counterclockwise from s — 
m 2 + iO to s — m 2 — iO around the cut of the perturbative 
polarization function n pt (including unphysical part). 

Siince 0^ +A is small, we shall be concerned with the 
D = 4, 6 operators. They can be conveniently separated by 
taking the real part of the Borel transform ( plj ) with com- 
plex argument: 
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At (j> = it/ '6 the operator 6 disappears and at <f> = 7r/4 
there is no O4 in the r.h.s. Fig. |]a,b displays both these pos- 
sibilities. Vertical bars correspond to Re B cxp , while the 
solid line with shaded area around it, labeled with "0.355" 
mark, shows purely perturbative contribution Re B pt com- 
puted with initial condition a s (m T ) = 0.355. The same 
is done for a s {m 2 ) = 0.330, which is the lowest possible 
value within the error range ( |T^ ) (for a reason which will 
become clear later). 

Consider at first <fi = ir/4. The gluonic condensate 
(G 2 ), and consequently O v+A must be positive. There- 
(19) fore, the perturbative curve must go below experimental 
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Figure 5: Sum rule ( p2[ ) with <p — (a) and (f> = 7r/4 (b). The dash line is the contribution of the gluonic condensate 
equal to the central value of (E3) added to the 0.330-perturbative curve. 



one, if the discrepancy is explained by OPE. However 
Fig. ||a shows, that the central value of 0.355-theoretical 
prediction goes above the experimental band for M 2 < 
0.9 GeV 2 . If we forget for a moment about theoretical 
uncertainty and assume that the theory should work for 
M 2 > 0.6 GeV 2 , as follows from our V — A analysis, 
this means that the condensate must be negative. It rather 
contradicts to our expectations. 

But if we take slightly lower input a s (to 2 ), the perturba- 
tive curve will go down. The lowest possible value is 0.330, 
as follows from (|l8]). Indeed, in this case the central theo- 
retical value is below experimental bars. If one takes some 
point, say, M 2 = 0.8 GeV 2 , then the following values of 
the gluonic condensate are acceptable: 



^ G »* G %) = U.i)0(ixf).0J2GcV' 



for a s (m 2 T ) = 0.330 and Af 2 >0.8GeV 2 (23) 

Theoretical and experimental errors are added together. In 
principle our result ( |23| ) does not contradict SVZ value. 
However, in order to achieve it, we must sit at "the very 
edge of errors", which seems unlikely. 



is posi- 



Now we turn to = tt/4, Fig. ]5|b. If 0\ +A 
tive, as OPE+factorization predict (EQ), its contribution to 
the r.h.s. of ( p2[ ) must be negative. This however strongly 
disfavors our previous choice a s (m 2 ) — 0.330, motivated 
by the sign of the gluonic condensate. So it seems rather 
difficult to make both (G 2 ) and D — 6 operator positive 
simultaneously. Here the OPE predictions are in certain 
disbalance with the data. 

This is not however a serious disagreement. First, both 
OX +A and 0^ +A are small enough. Available theoretical 
and experimental accuracy is about 2 — 3%, which is not 
sufficient to specify the values of both operators or to say 
something definite about their signs. 



Second, these operators are not rigorously defined ob- 
jects in perturbation theory. We do not have an algorithm 
to find their values from the first principles. So any state- 
ment about their properties should be considered with care. 
Moreover, it is not clear whether one could define them in- 
dependently of the perturbative series, to which they are 
added. In fact, different prescriptions to separate the so- 
called perturbative and nonperturbative terms may lead to 
different results, at least at the level of 1%. 



5 SUMMARY 



1. The V — A polarization operator Iiy\ A {Q 2 ) is well 
described by OPE series for Q 2 > 1 GeV 2 . 

2. The operator 0^~ A is approximately 2 times larger 
than what expected from QCD and low energy theo- 
rems. 



3. The V 



A polarization operator Ily V2 (Q 2 ) is well 



described by purely perturbative terms for Q 
lGeV 2 . 



> 



4. Current theoretical and experimental accuracy is not 
sufficient to determine the value of the gluonic con- 
densate (G 2 ). However it is likely to be much lower 
than commonly accepted SVZ value. 
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